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Motivated by the recent experimental realization of a candidate to the Fulde-Ferrell (FF) and 
the Larkin-Ovchinnikov (LO) states in one dimensional (ID) atomic Fermi gases, we study the 
quantum phase transitions in these enigmatic, finite momentum-paired superfluids. We focus on 
the FF state and investigate the effects of the induced interaction on the stability of the FFLO 
phase in homogeneous spin-imbalanced quasi-lD Fermi gases at zero temperature. When this is 
taken into account we find a direct transition from the fully polarized to the FFLO state. Also, we 
consider the effect of a finite lifetime of the quasi-particles states in the normal-superfluid instability. 
In the limit of long lifetimes, the lifetime effect is irrelevant and the transition is directly from the 
fully polarized to the FFLO state. We show, however, that for sufficiently short lifetimes there is 
a quantum critical point (QCP), at a finite value of the mismatch of the Fermi wave-vectors of the 
different quasi-particles, that we fully characterize. In this case the transition is from the FFLO 
phase to a normal partially polarized state with increasing mismatch. 
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I. INTRODUCTION 



Recently, there has been increased interest in the theory of one-dimensional (ID) imbalanced Fermi systems, partly 
because of the relevance of these theories for the understanding of the Fulde and Ferrell [l[ and Larkin and Ovchin- 
nikov @ (FFLO) phase. The FFLO is an exotic phase proposed approximately forty years ago, where atoms of opposite 
momenta and spins form Cooper pairs with finite momentum. In spite of intense theoretical and experimental efforts, 
the FFLO phase remains elusive. 

In three-dimensional (3D) systems, in the strongly-interacting limit, experiments show [H-Q that the gas phase 
separates with an unpolarized superfluid core surrounded by a polarized shell [ill, [HI, with no evidence for the 
FFLO phase 0. However, in one-dimensional (ID) imbalanced Fermi gases, the observed density profiles [8( agree 
quantitatively well with theories that exhibit the ID equivalent of FFLO correlations at low temperatures [9l4l2|. 
These experimental measurements @ of density profiles of a two spin mixture of ultracold ®Li atoms trapped in 
an array of ID tubes show that at finite spin imbalance, the system phase separates with an inverted phase profile 
as compared to the 3D case. In these ID experiments a partially polarized core was observed surrounded by wings 
composed of either a completely paired or a fully polarized Fermi gas, depending on the degree of polarization. 

This recent experimental observation of what can be seen as a strong candidate for FFLO-like correlations in ID 
has motivated theoretical investigations of possible mechanisms responsible for its stability. The increased stability of 
FFLO-like phases in ID can be understood as a nesting effect, where a single wavevector connects all points on the 
Fermi surface, allowing all atoms on the Fermi surface to participate in finite momentum pairing, while in 3D, only a 
small portion of these atoms are able to contribute. This ID Fermi surface nesting enhancement of the instability of 
the normal to the FFLO state is analogous to the conventional charge density wave (CDW) instability (l5| . 

In the paring mechanism, besides the particle- particle channel considered by Nozieres and Schmitt-Rink (NSR) |l6j , 
there is a correction of the two-body pairing interactions considered first by Gorkov and Mclik-Barkhudarov 
(GMB) [17j . This correction accounts for the induced interactions which arise between atoms at the Fermi level 
due to the polarization of the medium. It has been shown that these induced interactions suppress the superfluid 
transition temperature by a factor of about 2.22 in 3D [l8| and 2.72 in 2D [2(J spin-balanced Fermi gases, 
respectively, when compared with the mean-field (MF) results. The GMB correction was considered recently in 
various situations as, for instance, in a spin-balanced Fermi gas in an optical lattice [2lL [2~2| . in a homogeneous 
three-components Fermi gas [23[, and in the unitary limit of spin-balanced [24[ and imbalanced 3D Fermi gases (25j . 

In this paper we study the zero temperature (T) phase diagram of a quasi-ID imbalanced Fermi system as a function 
of the mismatch h between their Fermi wavevectors. This is relevant for the nearly ID Fermi gases we are interested 
in. We assume the existence of weak inter-tube interactions so that the MF BCS approach is justified at zero T (26j . 
We show that including induced interactions is essential to correct the MF results, since they give rise to a finite 
critical field (or mismatch) separating a fully polarized phase from the FFLO state, and in this way, to reveal the 
presence of the nesting effect in spin-imbalanced quasi-ID (ideal) Fermi gases. 

It is possible to conceive in actual physical systems mechanisms by which the quasiparticlc states in the normal phase 
acquire a finite lifetime. For example, due to different types of unknown (or unrecognized) scattering mechanisms not 
included in the pairing interaction or artificially engineered disorder. We show here how this lifetime effect modifies 
the T = phase diagram of the ID gas. In the limit of short lifetimes of the quasiparticle states in ID we find a 
quantum phase transition from the normal-to-inhomogeneous superfluid as the Fermi wave- vector mismatch is reduced 
from the normal phase. This T = transition is continuous or second order. It is associated with a QCP at a critical 
value of the field (mismatch) h c . We fully characterize this QCP obtaining its dynamic quantum critical exponent 
and universality class. On the other hand, for sufficiently long lifetimes (weak disorder), lifetime effects turn out to 
be irrelevant and we recover the previous results of including only induced interactions. Our results imply that for 
sufficiently strong disorder, the region in the phase diagram where the FFLO phase exits is reduced without necessarily 
being destroyed, even in ID. The instability of the normal state that we consider is that for a FFLO superfluid state 
characterized by a single wave-vector q. This is the first zero temperature instability that occurs as the effective 
Zeeman field h is reduced (26j from the normal phase. 

The majority of the recent literature on ID spin-imbalanced Fermi gases considers the simplest possible system that 
exhibits FFLO-type pairing namely, the Yang-Gaudin model or its lattice version, the Hubbard model with attractive 
interactions (27H30| . Both models are exactly solvable (or intcgrablc) and their energy spectra and thermodynamical 
properties can be calculated exactly using the Bethe ansatz and numerical methods [23]. However, although the Bethe 
ansatz yields exact solutions there is no consensus whether this approach can completely elucidate the nature of the 
ID FFLO state 13111 . Also out of reach of Yang-Gaudin model are inter-tube coup ling s so that one has to engage MF 
approximation [29( and beyond, as in the present work, or effective field theories 32 [. 

The theoretical predictions using the integrable Yang-Gaudin model [HI [H| agree very well with the experimentally 
observed density profiles and support the stability of the FFLO phase in ID. Here we give a robust and transparent 
physical explanation of this stability as a consequence of nesting effects. Wc also show explicitly that the medium 
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indubitably modifies the fermion-fermion interaction g, due to many-body effects. This effect is not captured by the 
Yang-Gaudin model. We point out that our results including induced interactions are consistent with the exact results 
since, for long lifetimes, we find a direct transition from the fully polarized to the FFLO state. 

Our study is based on calculations of pair and density fluctuations, i.e., of the particle-particle (pair) and particle- 
hole susceptibilities, respectively, that certainly are present in a Fermi gas, in any dimension. In ID these two physical 
quantities are of extreme importance, since both diverge, which are indications of some order in the system. As will be 
shown below, our main results do not depend on the particular Hamiltonian used to describe the attractive fcrmionic 
gas. Rather they arise from fundamental quantities, namely the ID fermionic dispersion relations that always can be 
linearized close to the Fermi points, regardless of its precise nature, and the subsequent calculation of the particle- 
hole and particle-particle susceptibilities [33} • Thus, given that it is of fundamental importance to consider the Fermi 
surface properties at ID or, the particle-particle and particle-hole interactions near it, our results can be considered 
as model independent. In this sense, our approach is complementary to those previous studies based on the real space 
Yang-Gaudin and Hubbard models. 

II. MODEL HAMILTONIAN 

To begin, let us consider a non-relativistic dilute (i.e., the particles interact through a short-range attractive 
interaction) ID spin-polarized Fermi gas, described by the following single-channel model Hamiltonian 

H = H-*^2n a n a (1) 

k . a 

= £ k a l a k + 4 b k b k + 9lD a k' b -k' b -k a k, 
k k,k' 

where at, etfc are the creation and annihilation operators for the a particles (and the same for the b particles) and ej* 
are their dispersion relation, defined by ej* = — pb a , with = h 2 k 2 /2m and [i a the chemical potential of the non- 
interacting a-species, a = a, b. A special case described by this Hamiltonian is that of identical spin 5 = 1/2 particles 
under an external magnetic field h. In this case, the a and b correspond to the spin up and spin down bands with 
their degeneracy raised by the magnetic field. The dispersion relations are then given by e^' b = — [i a ,b — £fc — M~F h. 
In particular, we will approximate the dispersions around their Fermi energies by e^ b = vp(k — kp) =F h, since the 
relevant states are near the Fermi momenta, vf is the Fermi velocity. For future notation, we rewrite this equation 
as e^' b = vp(k — k F b ), where k F ' b = kp ± h/vp are the Fermi wave-vectors of the different particles or bands. An 
important parameter in the present study is the Fermi wave-vector mismatch, Skp = k F — k F = 2h/vp. To reflect an 
attractive (s-wave) interaction between particles a and b we take gw = 9 < 0. 

III. THE IDEAL CASE: INFINITE LIFETIMES 

In this section we consider quasi-lD imbalanced Fermi system as a function of the mismatch h between their 
chemical potentials in the ideal case, i.e., in the absence of lifetime effects. Our aim is to investigate the quantum 
phase transition from the normal-to-inhomogeneous FFLO superfluid phase as the Fermi wave-vector mismatch is 
reduced from the normal polarized state. 

A. Ginzburg-Landau Theory and the FFLO Phase 

Since in a homogeneous ID system the quantum pha se transition to the FFLO phase is continuous [ll|, EH EBI j we 
expand the action in fluctuations, |A^|, and obtain |34l. |35| . 

Saff = S)/ ^o(a(|« ? 1, Wo )|A ? ( Wo )| 2 +0(|A| 4 )), (2) 

where cc(|g|, ojq) — — -— Xpp(|^|, Wo); with (q, ojq) being the external momemtum of the particle-particle bubble diagram 
and Xpp(|<7|, wo) the pair susceptibility, 
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XpP (\q\,u ) = cb q ' — 
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(3) 



(4) 



£<; + /i + (a + wq/2) a; + /i — (a + ujq/2) 



where cj c is an energy cut-off and a = vp\tf\/2. In the zero temperature limit, XppGsIj w o) is given by 
Re X (q,w ) =iV(0) 



In(-^j--ln|l-(? + £Dor|-jln|l-(?-wo) 



(5) 



where iV(0) = is the density of states for both spins at the Fermi energy Ep = k F /2m, q = vp\q\/2h, ujq = u>o/2h. 
Here vf is the Fermi velocity, and 



Im x(q,wo) 



-iV(0)|, if ft < (a - lo /2)\ 

- N(0)^, if (a - uj /2) 2 <h<(a + uj /2) 2 , 

, iih> (a + wo/2) 2 . 



(G) 



The static pair susceptibility reads 



XpP {q) = N(0) 

At the continuous phase transition, we find a{q) = 



-- — Xp P (<l) = 0, which yields the critical field h c (q): 



(7) 



fee 

Ao 



2^1 



(8) 



where Aq = 2w c exp(— 1/N(0)\g\) is the zero temperature BCS gap. This expression diverges for q = q c = 1 i.e., for 



— , yielding h c = oo. Since the two species Fermi momenta can be written as k F ,a = kp ± — , their difference 



k F = and we find that at the critical mismatch the pair wave vector reads q c 



k F . Thus, the 



q c : 
is k b F 

calculation above shows that FFLO type of correlations are so strong in quasi-lD systems with attractive interactions, 
that differently from d = 2 [37} and d = 3 [38|, in ID the FFLO phase persists even in the presence of an arbitrarily 
large magnetic field [36J. However, we show below that including induced interactions substantially modifies this 
scenario and there is a finite hf beyond which FFLO correlations disappear. This field coincides with that at which 
the system becomes fully polarized in agreement with the exact results [32j. 



B. Induced Interaction in a Spin Polarized Fermi Gas 



NSR have shown that as the superconducting transition temperature T c is approached from above, Cooper pair 
fluctuations grow in amplitude, and the pair susceptibility (which measures the tendency of pairs to form in response 
to an external pair field) diverges. The pair fluctuation is expressed by the pair susceptibility of Eq. (j3]),and NSR 
showed that at(\q\ = 0, u>o = 0) = is simply the Thouless condition for weak coupling superconductivity [16j|. 

Besides the pairing fluctuations that must be taken into account, as pointed out by NSR, to obtain the correct 
superfluid transition temperature T c of the BEC-BCS crossover, there is another effect of particle-hole fluctuations 
that affects the superfluid state. Namely, there is a change in the coefficient a(|g|) due to screening of the interspecies 
(or induced) interaction, known as the GMB correction [17|. In the BEC side, the NSR fluctuation is dominant, while 
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FIG. 1. The lowest order particle-hole diagram which generates the induced interaction <7md(pi,P4)- Arrowed and dashed lines 
describe fermionic propagators and the coupling g between the atoms, respectively. 

the GMB fluctuation becomes weaker towards the BEC side and vanishes in this region due to the disappearance of 
the Fermi surface (20j . 

In the original work by GMB [l7|, the induced interaction was obtained in the BCS limit by second-order pertur- 
bation (T7L EH . The diagram in Fig. [U describes a scattering process in which the conservation of total momentum 
implies that pi + P2 — > P3 + Pi- Besides, the frequencies and momentum are set to zero, so that p\ = —p2 and 
p*3 = — p^ This leads to the induced interaction giud{pi,Pi) = —g 2 Xph(pi — Pi) 0, where Xph(pi — Pi) is the 
particle-hole susceptibility, pi = (p^w^) is a vector in the space of wave- vector p and fermion Matsubara frequency 
lui = (21 + l)n/(hf3), p — l/(/cgT). The term gind{pi,Pi) describes the modification of the inter-particle interaction g 
due to many-body effects (i.e., the presence of other particles). 

In general, if one is dealing with small interaction g or the number of electrons in the band is such that the particle- 
hole susceptibility is well-behaved, the second order term in the induced interaction is sufficient to describe the system 
and calculate the change of the critical parameters. However, for general cases the GMB treatment can be extended 
to a region where all orders are important, and one has to sum an infinite series of the diagram represented in Fig. [T] 
in a random phase approximation (RPA). This yields [24l. |25|: 



g-md(pi,Pi) 



g 2 x P h(pi 



Pi) 



i + gx P h(pi - Pi)' 

Including the total induced interaction, the effective interaction between atoms with different spins is given by: 



(9) 



g e s(pi,Pi) = Syr = g + g-md(pi,Pi) = 



1 + gx P h(pi,Pi) 



(10) 



Next, we calculate the polarization function Xph(pi — Pi) = Xph{p') for a ID spin polarized fermionic gas. This is 
given by, 



X P h(p') = ■^^2Gv,b(jp)GaAP + P') = J 



dk 



flab — /k+k'a 



(2tt) 2 iH^ + e b k - e£ +k , 



(11) 



where p' = (k',Qi), fi; = 21tt/(H(3) is the Matsubara frequency of a boson, L is the size of the system, f(k) is the 
Fermi distribution function f(£^' b ) = l/(e /3£ '= + 1), with (3 = 1/T, where we have set ks = h = 1- The Matsubara 
Green's function of a non-interacting Fermi gas is given by Go<j(p) = l/(*w/ — e^ a ). We calculate the polarization 
bubble in the static limit (f2; = 0) and find 



X P h(k') = -N(0)f(x,h), 



(12) 



where k' = \k'\ is equal to the magnitude of p± + = p\ — so k' = y/(pi +P3)-(pi + P3) = \J p{ + P3 + 2pi .p>3 = 
\Jp\ +P3 + 2|pi||p3| cos</>, where <f) is the angle between pi and P3. Since the scattering is in ID, the only values of 
</> are or 7r. Performing the calculation we obtain the real function /(x, h) in the form of a generalized Lindhard 
function given by, 



In 



1 + x b H b 



rj b H b 



In 



1 + x a H° 



1 - x a H a 



(13) 
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(a) 



k a k k b 

F F K F 



(b) • 

b 

k F,eff 

FIG. 2. Geometrical illustration of the Fermi surfaces of a ID imbalanced gas for a given mismatch 5k F = k h F — kp. In (a) we 
see the Fermi points and wavevectors of the a and b species, relative to kF, and the nesting vector k! = kp + kp = 2&f. In (b), 
the system is fully polarized and non- interacting, due to the induced interaction, with an effective b- Fermi surface kp^ff = 2&f, 
with kp — 0. This fully polarized state occurs for a field hf = fi = ^uf^f. 

where x = 2^:, x b ' a = -X^, k F a = y/2mfi b ' a , H b ' a = l±ffi. Note that at h = 0, the well known result for ID 
balanced systems is recovered, 



/(x) = ^m 



(14) 



for which Xph{k' = 0) = —N(0). The function f(x,h) diverges for both x b ' a H b ' a = 1 or for either, x a H a = 1 or 
x b H b = 1. Solving, for example, x b H b = 1 for k' we find: 



k' h = k a F + k b F = 2k F . (15) 



Thus, at k' = k' = k F + k F = 2k F the function f(x, h) in Eq. (fT5)) diverges (a similarly condition holds for k' a ). This 
corresponds to = and = \ps\ = kp = \/2M[i, meaning that both scattering particles are at the Fermi surface. 
Notice that in metallic systems, the divergence of Xph(k') for a given value of g is often related to an instability to 
a charge ordered phase [2l[. The result above shows that even in the presence of an external magnetic field, the 
particle-hole susceptibility diverges at the same value of k, as in the absence of the "field" h, namely for k' = 2kp. 

Considering particle-hole fluctuations, we replace g by g c ff and the instability condition for the superconducting 
phase is given by: 

a eff(q) = -l - 9eff(k')Xpp(q) = o (16) 

It can be easily verified from this equation that the wave vector for which this condition is first satisfied is still given 
by Q = Q c = lj which gives a critical field h c = oo. An interesting and new possibility occurs when kl = k b F + k F = 
2kp = q c = k F — k F . However, this is possible if, and only if, k F = 0. Then, we conclude that the many-body effects 
brought about by the nesting wave vector k' which connects the two Fermi points k F and k b F gives rise a new effective 
b species Fermi surface, k b F e ^ = 2kp, with k F = 0. Indeed, it can be seen from the divergence of the particle-hole 
susceptibility and Eq. (|10[) that the effective interaction is <? e // = 0. This situation q c = 2kp, as we just verified, 
corresponds to the fully polarized gas. This fully polarized gas, which is equivalent to empty the band of down spins 
and accommodate all in band 6, such that, k F = and k F = 2fcp, is reached for a field hf = /x = ^vpkp, as illustrated 
in Fig. [21 Then for h > hf the fully polarized system can be considered as non-interacting and it remains normal for 
h > hf. For h < hf the system enters the FFLO phase. 

b _ a , 

The spin polarization is defined as P = ™b + ™ a > where n a ' are the number densities. Since in a ID system we have 
kp b = ^n a ' b , the polarization can be written as: 

p _ (^ + M 1/2 -(^-M 1/2 r . 

(ji + hyP + Qi-h) 1 /*' [ ' 

For hf = fx, we have P = 1 as expected. 

We conclude this section with the result that in the ideal case, where the quasi-particles have an infinite lifetime, 
the FFLO phase will occur for all h < hf and for any strength of the attractive interaction. Here hf is the field above 
which the system is fully polarized. The quantum phase transition in this case is directly from the fully polarized 
state to a phase with FFLO correlations (321 ] - 
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FIG. 3. (Color online) Critical magnetic field for the normal-to-helicoidal superconductor quantum phase transition as a function 
of normalized wave-vector for pure (full line) and homogeneous, with finite lifetime, (dashed line) quasi-one dimensional systems. 



IV. LIFETIME EFFECTS 



In this section we study lifetime effects in the phase diagram of ID attractive Fermi gases. These arise whenever 
the wavevector k is not a good quantum number to describe the quasi-particles states. For example, it is not difficult 
to imagine engineering disorder in an optical lattice, to break translational invariance, such that the quasi-particle 
states acquire a finite lifetime [39^ . 

In this case it is appropriate to introduce an imaginary self-energy in the dispersion relation of the fermions. Taking 
into account the finite lifetime of the quasi-particle states, the particle-particle dynamic susceptibility can be written 
as, M 



N(0) 



dx 



1 



+ 



x + h + v F q/2 - uj q /2 + ij/2 x + h- v F q/2 + uj q /2 - ij/2 



+ (v F q/2 -> -v F q/2) 



(18) 

where 7 = r 1 , the inverse of the lifetime of a quasiparticle state in the normal phase. This approach is formally 
similar to that used to investigate the effect of non-magnetic impurities in higher dimensional (3D and 2D) FFLO 
superconductors [ill |42| . In this case the main interest was to obtain the reduction in the critical temperature of the 
superconductor. Here we will concentrate in the zero temperature phase diagram of the ID system and will be able 
to fully characterize the new QCP that arises due to the lifetime effects. 

Since we are interested in studying the effect of this finite lifetime on quantum criticality, we start calculating the 
real part of the static particle-particle susceptibility 



Kex(g, w = 0) 



N(0) 



dx 



x + h + v F q/2 



+ h — v F q/2 



(x + h + v F q/2) 2 + 7 2 /4 (x + h - v F q/2) 2 + 7 2 /4 



This can be easily integrated and for uj c t >> 1, we get, 

Vlex(q, loo = 0) = N(0) In^-^P- In 

h 4 



The condition for the divergence of the interacting pair susceptibility becomes, 



ffr^(0) 



2h 1 

in h — In 

A 4 



(l + q) 2 + 



21 r 



7 
4/i 2 



(W) 2 + 



2 1 



J 

Ah 2 



0, 



(19) 



(20) 



(21) 



where g r is an effective (renormalized by the GMB correction) coupling constant. This condition is first satisfied when 

maxii 

2h{q) 



the argument of the logarithmic is maximum, that is, for q c = j 2 /4h 2 \. The critical field is given by h c = h(q c ) 

where. 



1 



An 



:i + g) 2 + 3& (W) 2 + i& 



1/4' 



(22) 



Substituting for q c we finally get, 



h(q c ) = Ag/47, 



(23) 
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implying that in quasi-lD imbalanced Fermi systems, for sufficiently short lifetimes of the quasiparticle states, the 
FFLO phase appears only below a critical mismatch h < h c = Aq/4-j. 

Taking into account the results of the previous section we can summarize our results as: 

• If h c > hf, where hf was calculated in the previous section, the lifetime effect is irrelevant and the phase 
transition occurs at hf directly from the fully polarized state to the FFLO phase. 

• For short lifetimes, or strong disorder, such that, h c < hf, the transition to the FFLO phase occurs from a 
partially polarized state. In this case there is an intervening normal partially polarized (PP) phase between the 
fully polarized (FP) state and the FFLO phase (h c < h < hf). 

We emphasize that the results above do not arise from a MF calculation, but rather on the properties of the pair 
and density fluctuations that go beyond the MF approximation. Besides, they are not dependent on the Hamiltonian 
given byEq.(l), and rely on the ID fermionic dispersion relations that can always be linearized close to the Fermi 
points [33f . 



A. Nature of the transition at h c < h 



f 



From here on, we will consider the case h c < hf and study the quantum phase transition from the normal PP state 
to the FFLO phase. The real part of the dynamic susceptibility is obtained as, 



N(0) 



In 



)2 + i^ 

2 



(1 + + i L- 



(1-9- loo) 2 
(1 + g + oJo) 2 - 



7 
" 4h 2 

Ah 2 



(24) 



Expanding close to q — q c and uq = 0, we get 

l-g r Viex(q,Lu )~ g r N{0) 



In- 



h c ' 7 

The imaginary part of the dynamic susceptibility is given by, 



+ (^) 2 (g-g c ) 2 + (^) 2 cJ 2 



(25) 



3mx(<?,w ) = 



N(0) 



dx 



7/2 



7/2 



[x - (w /2 - h - v F q/2)f + 7 2 /4 [x - (wo/2 -h + v F q/2)] 2 + 7 2 /4 



+ (luq -> -w ) 



which can be easily integrated to give, 
^mX{q,u} ) = — 



tan 



2(q-u ) 



tan 



2(q + uJ ) 



7 1- 



7 1- 



(q+LJoY 2 — 1 



(26) 



(27) 



where we recall, q = VFq/2h, ujq = u)o/2h and we defined 7 = 7/2/1. 
Expanding close to uj = 0, q = q c we obtain: 



Smx(<7,w ) = N(0)-u} . 

7 



(28) 



Then in the limit ujq — > 0, the frequency dependent part of the imaginary susceptibility dominates over the real part. 

As concerns its quantum critical behavior, the zero temperature phase transition from the normal to the inhomoge- 
neous superconductor state of the homogeneous quasi- ID system in the presence of a finite lifetime of the quasi-particle 
states can be described by the following effective action, at the Gaussian level [Ifl , 



S eff = JdQ J duj [S + Q 2 + \uj a \] |A(Q,wo) 



(29) 



where 5 = h — h c , and Q = q — q c . The quantum critical point associated with this phase transition has a dynamic 
exponent z = 2, such that, its effective dimension, d e ff = d + z = 3 (43l |44|. Consequently, we expect the superfluid 
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transition in the ID system in the presence of a finite lifetime to be in the universality class of the 3D XY model [44 1 
due to the two-component nature of the superfluid order parameter. Effects of temperature [45j can also be obtained 
from knowledge of the critical exponents of the quantum critical point. In this case the finite temperature critical 
line is obtained as T c oc \h — h c \ vz where v rj 2/3 is the correlation length exponent of the 3d-XY model (4(| and 
z = 2 as obtained previously. Since there is no long range magnetic order in d = 1 at finite temperatures, this line in 
practice provides the temperature scale below which the FFLO correlations become important and this varies with 
the distance to the QCP. Since (dT c /dh)h=h c = 0, the characteristic temperature T c turns out to be very small near 
the critical field. 

The identification of the universality class of the lifetime induced QCP as being that of the 3d XY model also allows 
to obtain the behavior of the correlation function of the FFLO fluctuations at T = 0. Usin g th at the exponent 77 for 
the order parameter correlation function of the 3d XY model takes the value rj = 0.0381 [43], we find that at the 
QCP the FFLO correlation function decays with distance r, as G(r) cx i/ r d + z ~' 2 +n [44^ i.e., G(r) = i/ r 1 038 . This 
exponent turns out to be small or of the same order of that obtained numerically for these type of correlations using 
the Bethe-ansatz (see Refs. [28[ and [32l]). 



V. DISCUSSION AND CONCLUSION 



In spite of theoretical predictions and intensive experimental activity, the FFLO phase remains elusive. Motivated 
by experimental results in cold atom systems and aiming to understand the reasons for the difficulty in observing 
this phase, we have carried out a detailed study of a FFLO phase in ID systems, which provide the most favorable 
conditions for the appearance of this phenomenon. We have considered the effect of the induced interaction in an 
ideal gas to show how the mean-field result which predicts long range FFLO correlations in the presence of arbitrarily 
large magnetic field or mismatches is modified leading to a finite critical field hf . At this field, a nesting condition is 
satisfied and long range FFLO correlations set in. In agreement with the exact results, the system at hf goes directly 
from a regime with strong FFLO correlations (h < hf) to the fully polarized normal phase (h > hf). 

In systems, with additional interactions not included in the pairing Hamiltonian, or in the presence of artificial 
disorder, a finite lifetime of the quasi-particle excitations may be considered. We have shown here how the lifetime 
effect modifies the T = phase diagram in ID. It gives rise to a new characteristic or critical field h c which depends 
on the h = BCS gap and on the lifetime, t = I/7 of the states (Eq. [23]) . If disorder is weak, such that, h c > hf, 
disorder is irrelevant and the transition with increasing field is from the FFLO phase to the fully polarized system at 
hf. However, for strong disorder (h c < hf) there is a new QCP in the system that we have fully characterized. In this 
case with increasing mismatch the system goes from the FFLO phase to a normal partially polarized phase at h = h c 
and finally to a fully polarized phase at hf. We have fully characterized the QCP at h c . In this case the region of the 
phase diagram where the FFLO phase appears is reduced. 

We hope our results will stimulate further experiments to confirm unambiguously the existence of the FFLO phase 
in quasi- ID imbalanced Fermi gases. 
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